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Abstract 


To trace the evolution of canine coronavirus (CCoV), 201 stool samples from diarrheic dogs 
in northeast China were subjected to reverse transcription-polymerase chain reactions (RT- 
PCRs) targeting the partial M and S genes of CCoV, followed by an epidemiological analy- 
sis. M gene RT-PCRs showed that 28.36% (57/201) of the samples were positive for CCoV; 
of the 57 positive samples, CCoV-I and CCoV-II accounted for 15.79% (9/57) and 84.21% 
(48/57), respectively. A sequence comparison of the partial M gene revealed nucleotide 
homologies of 88.4%—100% among the 57 CCoV strains, and 88.7%—96.2% identity 
between the 57 CCoV strains and the Chinese reference strain HF3. The CCoV-I and 
CCoV-II strains exhibited genetic diversity when compared with reference strains from 
China and other countries. The 57 CCoV strains exhibited high co-infection rates with 
canine kobuvirus (CaKV) (33.33%) and canine parvovirus-2 (CPV-2) (31.58%). The CCoV 
prevalence in diarrheic dogs differed significantly with immunization status, regions, sea- 
sons, and ages. Moreover, 28 S genes were amplified from the 57 CCoV-positive samples, 
including 26 CCoV-lla strains, one CCoV-IIb strain, and one CCoV-I strain. A sequence 
comparison of the partial S gene revealed 86.3% —100% nucleotide identity among the 26 
CCoV-lla strains, and 89.6%—92.2% identity between the 26 CCoV-lla strains and the Chi- 
nese reference strain V1. The 26 CCoV-Ila strains showed genetic diversity when com- 
pared with reference strains from China and other countries. Our data provide evidence that 
CCoV-I, CCoV-Ila, and CCoV-IIb strains co-circulate in the diarrhoetic dogs in northeast 
China, high co-infection rates with CaKV and CPV-2 were observed, and the CCoV-II 
strains exhibited high prevalence and genetic diversity. 
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Introduction 


Canine coronavirus (CCoV) was first recognized as an enteric pathogen of dogs in 1971 [1]. 
CCoV is a common infection in young dogs, particularly those housed in large groups [2-5]. 
CCoV is an enveloped, single-stranded, positive-sense RNA virus, and it belongs to the family 
Coronaviridae, subfamily Coronavirinae, genus Alphacoronavirus, species Alphacoronavirus-1 
[6]. CCoV consists of two distinct genotypes, CCoV-I and CCoV-II; the CCoV-II viruses are 
further divided into two subtypes CCoV-I]a and CCoV-IIb [3, 7]. The S protein of CCoV is a 
glycoprotein peplomer on the viral surface, and it plays an important role in the induction of 
neutralizing antibodies, specific receptor binding, and cell membrane fusion. Accumulating 
reports attributed that the increase in the severity of CCoV infections in dogs and the emer- 
gence of CCoV variants to potential recombination events within the S gene, which occur 
when a host is co-infected with different CCoV types [8-10]. Therefore, CCoV has received 
much attention as an emerging cause of infectious disease in dogs [4, 11-16]. 

CCoV infection is a leading causes of diarrhea in dog population in China. Wang et al. 
(2006) reported that CCoV-II infections were very common in domestic dog, fox, and rac- 
coon-dog populations in China [17]. Ma et al. (2008) reported the molecular characterization 
of the 9.36-kb 3’ region of the CCoV 1-71 strain [18]. Gao et al. (2009) reported the isolation 
and identification of a CCoV strain from giant pandas in China [19]. However, information 
about the epidemiology of CCoVs in China is not available in the past five years. In the current 
study, we conducted a molecular epidemiologic investigation of CCoV in Heilongjiang prov- 
ince, northeast China. Moreover, the genetic evolution and co-infection of the identified CCoV 
strains were analyzed. Our aim was to provide insights into the epidemiology and genetic 
diversity of the CCoV strains circulating in northeast China. 


Material and Methods 
Ethics Statement 


The animal experiment, sampling, was approved by the Animal Care and Use Committee of 
the Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences, China. 
The sampling and data publication also were approved by animal’s owners. The field study did 
not involve endangered or protected species. No specific permissions were required for loca- 
tions of samples because the samples were collected from public areas or non-protection areas. 


Sampling 

A total of 201 fecal samples were collected in the form of rectal swabs of dogs with diarrhea 
from animal hospitals in the Harbin, Daqing, and Mudanjiang districts of Heilongjiang prov- 
ince in northeast China from May 2014 to April 2015, using 3.5-ml commercial virus sampling 
tubes (YOCON Biological Technology Co. Ltd., Beijing, China). For all samples, animal age, 
animal breed, animal gender, collection date, and vaccination were recorded, respectively. Of 
the 201 samples, 141 were collected in Harbin, 20 were collected in Daqing, and 40 were col- 
lected in Mudanjiang. All rectal swab samples were stored at —80°C, and they were also used 
for etiological investigations in our other studies [20, 21]. 


RNA extraction 


After 1 mL of fecal samples was centrifuged at 1,500 x g for 10 min at 4°C, the supernatant of 
each sample was transferred to a 1.5-ml Eppendorf tube. Viral RNA was extracted from each 
sample using the TIANamp Virus RNA Kit (Tiangen Biotech Co., Ltd., Beijing, China) accord- 
ing to the manufacturer’s instructions. The extracted RNA were stored at -80°C. 
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Detection and sequence analysis of CCoV 


Molecular detection of CCoV was conducted using reverse transcription-polymerase chain 
reaction (RT-PCR) targeting a 409-bp fragment of the M gene of CCoV that was described by 
Pratelli et al. (1999) [22]. The S gene fragments used for CCoV genotyping/subtyping, includ- 
ing a 346-bp fragment of CCoV-I, a 694-bp fragment of CCoV-Ila, and a 370-bp fragment of 
CCoV-IIb, were amplified using the primers EL1F/EL1R, S5F/S6R, and CEPol-1/TGSP-2, 
respectively [23, 24]. Briefly, first-strand cDNA was synthesized using random primers (six 
nucleotides) using Moloney murine leukemia virus (RNase H-) reverse transcriptase (Novo- 
protein Scientific Inc., Shanghai, China) according to the manufacturer’s instructions. In this 
study, the EmeraldAmp PCR Master Mix (2x Premix) (TaKaRa Biotechnology Co., Ltd., 
Dalian, China), and the Applied Biosystems GeneAmp PCR System 9700 thermal cycler 
(Thermo Fisher Scientific, Waltham, MA, USA) were used for PCR amplifications of all target 
genes. Other PCR amplification conditions were performed according to the protocols 
described by Pratelli et al. (2004) [23] and Erles and Brownlie (2009) [24]. The purified PCR 
products were directly subjected to Sanger sequencing, and each sample was sequenced three 
times. Sequence analysis was performed using the EditSeq program in the Lasergene DNAS- 
TAR™ version 5.06 software (DNASTAR Inc., Madison, WI, USA). All nucleotide sequences 
generated in our study were submitted to GenBank under accession numbers KT192642- 
KT192698 for the M gene of the 57 CCoV strains, KT222969-KT222994 for the S gene of the 
26 CCoV-Ila strains, KT222995 for the S gene of the CCoV-I strain, and KT222996 for the S 
gene of the CCoV-IIb strain. 


Phylogenetic analysis 


For the phylogenetic analysis, partial sequences of the M and S genes of CCoV strains, includ- 
ing CCoV-I, CCoV-IIa, and CCoV-IIb strains, were retrieved from GenBank. To construct the 
phylogenetic trees, a multiple alignment of all target sequences was performed using ClustalX 
program version 1.83. Furthermore, phylogenetic trees of all target sequences were generated 
from the ClustalX-generated alignments by MEGA6.06 software using the neighbor-joining 
method [25]. Neighbor-joining phylogenetic trees were built with the p-distance model, 1000 
bootstrap replicates, and, otherwise, the default parameters in MEGA 6. 


Screening for canine enteric pathogens 


The 57 CCoV-positive samples were screened for CPV-2, CaKV, canine astrovirus (CaAstV), 
canine norovirus (CNoV), canine bocavirus (CBoV), and Group A-Rotavirus (CRV-A) by 
either PCR or RT-PCR, followed by sequencing, according to previously described protocols 
[26-30]. All nucleotide sequences generated in our study were submitted to GenBank, and the 
accession numbers of the target genes of co-infected viruses are shown in S1 Table. 


Results 
Investigation and genotyping of CCoV 


The characteristics of the 57 CCoV-positive dogs, the genotyping of 57 CCoV strains, and the 
amino acid substitutions of M protein of 57 CCoV strains are shown in S1 Table, and a further 
analysis of the 57 CCoV-positive samples is shown in Table 1. Nucleotide sequences of the par- 
tial M gene of the 57 CCoV strains identified in our study were shown in Supporting Informa- 
tion (S1 Fig). Of the 201 samples, 57 samples (28.36%) were positive for CCoV following 
RT-PCR amplification of the partial M gene, and the CCoV-positive rates of the Harbin, Daq- 
ing, and Mudanjiang districts were 26.95%, 50%, and 22.5%, respectively (Table 1). Of the 57 
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Table 1. Further analysis of the CCoV-positive samples. 


Numbers of sample 
Positive rate for CCoV 
Genotyping of CCoV based on the partial M gene 
I 
ll 
Vaccination* 
Yes 
No 
Collection data 
Jan. to Mar. 
Apr. to Jun. 
Jul. to Sep. 
Oct. to Dec. 
Age 
O0< Age <1M 
1M< Age <2M 
2M< Age <3M 
3M< Age <4M 
Coinfection of CCoV with other enteric viruses 
Single infection 
CPV-2 
CakKV 
CBoV 
Mixed infection 
CPV-2+CakV 
CPV-2+CBoV+CakV 
Total co-infection 


CaKkKV, CPV-2, CBoV, CPV-2+CaKkV or CPV-2+CBoV+CakV 


Identity of partial M gene of CCoV identified in our study 
Nuleotides 
I 
Il 
I+l] 
Amino acids 
I 
Il 
I+ll 


Compared with M gene of the HF3 strain from China (AY864661) 


Nuleotides 
I+ll 
Amino acids 
I+ll 
Identity of partial S gene of CCoV identified in our study 
Nuleotides 
lla 
Amino acids 
lla 
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Harbin 


141 
26.95% (38/141) 


13.16% (5/38) 
86.84% (33/38) 


52.63% (20/38) 
34.21% (13/38) 


5.26% (2/38) 
2.63% (1/38) 
31.58% (12/38) 
60.53% (23/38) 


13.16% (5/38) 
42.11% (16/38) 
34.21% (13/38) 


36.84% (14/38) 
39.47% (15/38) 
2.63% (1/38) 


21.05% (8/38) 
2.63% (1/38) 


57.89% (22/38) 
95.7%-99.1% 
97.4%-100% 
88.4%-100% 
97.4%-100% 
100% 
94.8%-100% 


88.7%-96.2% 


92.2%-97.4% 


89.3%-100% 


92.1%-100% 
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Daqing 
20 
50% (10/20) 


10% (1/10) 
90% (9/10) 


20% (2/10) 
50% (5/10) 


60% (6/10) 


10% (1/10) 
30% (3/10) 


10% (1/10) 


50% (5/10) 
30% (3/10) 


20% (2/10) 
30% (3/10) 


50% (5/10) 


97.4%-100% 
92.5%-100% 


98.3%-100% 
93.9%-100% 


91.3%-95.9% 


93%-97.4% 


87.7%-100% 


87.6%-100% 


Mudanjiang 
40 
22.5% (9/40) 


33.33% (3/9) 
66.67% (6/9) 


33.33% (3/9) 
33.33% (3/9) 


22.22% (2/9) 
11.11% (1/9) 


66.67% (6/9) 


33.33% (3/9) 
22.22% (2/9) 
44.44% (4/9) 


22.22% (2/9) 
11.11% (1/9) 
22.22% (2/9) 


55.56% (5/9) 


97.4%-100% 
97.7%-100% 
89.6%-100% 


100% 
100% 
95.7%-100% 


89%-95.9% 


93%-97.4% 


Total 


201 
28.36% (57/201) 


15.79% (9/57) 
84.21% (48/57) 


43.86% (25/57) 
36.84% (21/57) 


17.54% (10/57) 
3.51% (2/57) 

22.81% (13/57) 
56.14% (32/57) 


1.75% (1/57) 
22.81% (13/57) 


36.84% (21/57) 
29.82% (17/57) 


31.58% (18/57) 
33.33% (19/57) 
5.26% (3/57) 


14.04% (8/57) 
1.75% (1/57) 


56.14% (32/57) 
94.8%-100% 
96.8%-100% 
88.4%-100% 
94.8%-100% 
98.3%-100% 
93%-100% 


88.7%-96.2% 


92.2%-97.4% 


86.3%-100% 


87.6%-100% 
(Continued) 
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Table 1. (Continued) 


Harbin Daqing Mudanjiang ___ Total 
Compared with S gene of the V1 strain from China (AY390342) 
Nuleotides 
lla 89.6%-92.2% 89.9%-92.2% — 89.6%-92.2% 
Amino acids 
lla 95.5%-97.5% 91.1%-97.5%  — 91.1%-97.5% 
Amino acid substitution of partial M protein of CCoV identified in our study 
'8Thr/Val 
| (Val—Thr) 100% (5/5) 0 (0/1) 100% (3/3) 88.89% (8/9) 
ll (Thr—Val) 100% (33/33) 100% (9/9) 100% (6/6) 100% (48/48) 
198 \Vlet/lle 
| (Ile—Met) 100% (5/5) 100% (1/1) 100% (3/3) 100% (9/9) 
Il (Met—lle) 100% (33/33) 100% (9/9) 100% (6/6) 100% (48/48) 
206His/Asn 
| (Asn—His) 100% (5/5) 100% (1/1) 100% (3/3) 100% (9/9) 
Il (His—Asn) 100% (33/33) 100% (9/9) 100% (6/6) 100% (48/48) 
228GIn/Lys 
| (Llys—Gln) 100% (5/5) 0 (0/1) 100% (3/3) 88.89% (8/9) 
Il (GIn—Lys) 100% (33/33) 100% (9/9) 100% (6/6) 100% (48/48) 


*The animals are vaccinated for CPV-2, canine distemper virus, canine parainfluenza virus, canine adenovirus type 1, and canine adenovirus type 2. 


doi:10.1371/journal.pone.0146975.t001 


CCoV-positive samples, CCoV-I and CCoV-II accounted for 15.79% (9/57) and 84.21% (48/ 
57), respectively; immunized animals and non-immunized animals accounted for 43.86% and 
36.84%, respectively; 56.14% (32/57) were collected from October to December, and 66.66% 
(38/57) were collected from dogs aged from 2 to 4 months; the total co-infection rate of CCoV 
with any viruses reached 56.14%, and CaKV and CPV-2 accounted for 33.33% (19/57) and 
31.58% (18/57), respectively (Table 1). Of the 57 CCoV strains, only 28 S gene sequences were 
successfully amplified, of which 26 belonged to CCoV-IIa, one to CCoV-IIb, and one to 
CCoV-I (S1 Table). Nucleotide sequences of the partial S gene of CCoV strains identified in 
our study were shown in Supporting Information (S2 Fig). 


Sequences and phylogenetic analysis of CCoV 


The sequence comparison of the partial M gene revealed nucleotide homologies of 88.4%- 
100% and amino acids homologies of 93%-100% among the 57 CCoV strains, while nucleotide 
and amino acid homologies of 88.7%-96.2% and 92.2%-97.4%, respectively, were observed 
between the 57 CCoV strains and the Chinese reference strain HF3. Of the 57 CCoV strains, 
the nine CCoV-I strains exhibited 94.8%-100% nucleotide identities and 94.8%-100% amino 
acid identities, and the 48 CCoV-II strains showed 96.8%-100% nucleotide identities and 
98.3%-100% amino acid identities (Table 1). A total of 12 amino acid substitutions were found 
in the partial M protein (S1 Table); the four substitutions, Thr178Val, Met198Ile, His206Asn, 
and Gln228Lys, occurred in all identified CCoV-II strains (48/48), and two substitutions, 
Ile198Met and Asn206His, occurred in all identified CCoV-I strains (9/9) (Table 1). The 
sequence comparison of the partial S gene revealed nucleotide homologies of 86.3%-100% and 
amino acids homologies of 87.6%-100% among the the 26 CCoV-IIa strains, and nucleotide 
and amino acid homologies of 89.6%-92.2% and 91.1%-97.5%, respectively, between the 26 
CCoV-IIa strains and the Chinese reference strain V1 (Table 1). 
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A phylogenetic analysis using the nucleotide sequences of the M gene revealed that the nine 
CCoV-I strains formed two clusters (C1 and C2), while the 48 CCoV-II strains formed six clus- 
ters (C1-C6). The 48 CCoV-II strains only exhibited a close relationship to one Chinese refer- 
ence strain, NJ17, and differed genetically from most of the reference strains from China and 
other countries (Fig 1). A phylogenetic analysis using partial S gene sequences demonstrated 
that the 26 CCoV-Ila strains were closely related to three reference strains, 5281 (Japan), TN- 
449 (USA), and 1086-IlIa (Brazil), and differed genetically from reference strains from China 
and other countries (Fig 2A). The CCoV-I strain HRB-A4 showed a close relationship to the 
reference strain from Italy (Fig 2B), while the CCoV-IIb strain HRB-BB9 was closely related to 
reference strains from European countries (Fig 2C). 


Discussion 


In our study, 57 of 201 samples (28.36%) were positive for CCoV, and the CCoV-positive rate 
differed among the three districts of Heilongjiang province (26.95% for Harbin, 50% for Daq- 
ing, and 22.5% for Mudanjiang). The CCoV-positive rate in feces has been reported to be 
43.75% (21/48), 32.83% (22/67), 50.46% (55/109), 54.39% (31/57), 12% (30/250), 51.3% (20/ 
39), 80% (8/10), 60% (3/5), 38.09% (8/21), 12.5% (1/8), 36.17% (17/47), and 11.54% (3/26) in 
China [17], the Republic of Korea [31], Japan [32], Albania [33], Brazil [11], Italy, Belgium, the 
Netherlands, Germany, the United Kingdom, Spain, and France, respectively [34]. These data 
demonstrated that CCoV infections showed clear differences among the geographical regions. 
The total CCoV-positive rate reported here was lower than that in a previous report in 2006 in 
China [17]. In our study, CCoV-II, accounting for 84.21% of the 57 CCoV strains, was the pre- 
dominant CCoV type in northeast China, which is line with most reports [11, 17, 34]. 

Mixed infections of canine enteric viruses frequently occur in diarrheic dogs. In our study, 
the total co-infection rate of CCoV with one or more CaKV, CPV-2, and CBoV strains was 
56.14%; single co-infections with CaKV, CPV-2, and CBoV accounted for 33.33%, 31.58%, and 
5.26%, respectively. Co-infections between CCoV and CPV-2 have been documented in dog 
populations in Western Europe, Japan, and Albania [32-34]. However, the high co-infection 
rate of CCoV with CaKV has not been reported in China and other countries. CaKV had been 
reported to be associated with severe enteritis in a litter of puppies [13]. It is speculated that the 
high prevalence of co-infections of CCoV with CaKV and CPV-2 may be associated with the 
occurrence of viral diarrhea in dogs in northeast China. Further studies should be conducted to 
understand the effect of the high co-infection rate with CaKV on the severity of clinical symp- 
toms of CCoV infections. 

In our study, CCoV-positive rate showed clear differences among seasons and ages. The 
high prevalence of CCoV in the feces of diarrheic dogs aged 2-4 months was partially validated 
in other studies [11, 32, 33]. The high positive rate (56.14%) of CCoV was found between Octo- 
ber and December, which may be associated with seasonal variations in northeast China. At 
present, most dogs in northeast China are vaccinated for CPV-2, canine distemper virus, 
canine parainfluenza virus, canine adenovirus type 1, and canine adenovirus type 2. In our 
study, of the 57 CCoV-positive samples, the immunized and non-immunized animals 
accounted for 43.86% and 36.84%, respectively. This result demonstrated that vaccination for 
other canine viruses did not effect CCoV infections in the dog population in Heilongjiang 
province, northeast China. Given the high prevalence and co-infection rates of CCoV, although 
it is controversial whether CCoV vaccines provide adequate immunity [35, 36], the immuniza- 
tion for CCoV is recommended in the dog population in northeast China in the future. 

Fragments of the M and S genes have been used to genotype CCoV [4, 11, 24, 37]. In our 
study, the CCoV-I and CCoV-II genotypes of the 57 CCoV strains were successfully identified 
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Fig 1. Phylogenetic analysis of CCoV strains using the nucleotide sequences of the partial M gene of 
CCoV. Red circles represent the CCoV strains identified in our study, and black diamonds represent the 
Chinese CCoV reference strains. 


doi:10.1371/journal.pone.0146975.g001 
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Fig 2. Phylogenetic analysis of CCoV strains on the basis of partial S gene sequences. (A) Phylogenetic tree generated using the nucleotide 
sequences of the partial S gene of CCoV-lla. (B) Phylogenetic tree generated using the nucleotide sequences of the partial S gene of CCoV-I. (C) 
Phylogenetic tree generated using the nucleotide sequences of the partial S gene of CCoV-IIb. Red circles represent the CCoV strains identified in our study. 
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using fragments of the M gene to construct a phylogenetic tree, coupled with an analysis of 
amino acids substitutions of the partial M protein. In the phylogenetic tree of the M gene, the 
48 CCoV-II strains identified in our study, which formed six clusters (C1-C6), differed geneti- 
cally from other country’s reference strains and the early Chinese strains (except strain NJ17), 
and they possessed characteristics of the CCoV strains in northeast China. Meanwhile, the 
nine CCoV-I strains identified in our study formed two clusters (C1 and C2). These data dem- 
onstrated that the CCoV-I and CCoV-II strains identified in our study exhibited genetic diver- 
sity. Additionally, of the nine CCoV-I strains, only one co-infection between CCoV-I and 
CCoV-II was found, accounting for 11.11% (1/9). This result was lower than that in a recent 
report in Brazil [11]. In the M protein, the Thr178Val, Met198Ile, His206Asn, and Gln228Lys 
substitutions were common in all identified CCoV-II strains when compared with amino acids 
substitutions in the partial M protein from all identified CCoV-I strains and reference strains. 
Co-infection of CCoV-I and CCoV-I strains has been documented in previous studies [17, 
32, 38-41]. To investigate single or multiple infections of CCoV types I, Ila, and IIb, we 
attempted to amplify three fragments of the S gene of 57 CCoV-positive samples by RT-PCR. 
Only the S genes of 28 CCoV strains were obtained, of which 26 strains belonged to CCoV-IIa, 
while the remaining two strains belonged to CCoV-I and CCoV-IIb, respectively. Compared 
with the reference strains, the partial S genes of the 26 Chinese CCoV-IIa strains exhibit genetic 
diversity, which is in agreement with the phylogenetic analysis of the M genes identified in our 
study. The CCoV-IIb strain HRB-BB9 was closely related to CCoV-IIb reference strains from 
Europe. The CCoV-I strain HRB-A4 clustered with reference strains from Europe, North 
America and South America. Additionally, it was rarely reported that one co-infection between 
CCoV-IIa and CCoV-IIb was found in the 28 CCoV-positive samples. The limited data 


8/11 


®-PLOS | ove 


Circulation of Canine Coronavirus in China 


demonstrated that CCoV-I, CCoV-Ila, and CCoV-IIb co-circulated in northeast China, and 
that CCoV-IIa strains predominated in the dog population. 


Conclusions 


The current study is the first to reveal that CCoV-I, CCoV-IIa, and CCoV-IIb strains co-circu- 
late in northeast China, and that there are high co-infection rate with CaKV and CPV-2. CCo- 
V-II strains are predominant, and they exhibit genetic diversity. Resulting data increase our 
understanding of the evolution of CCoV strains circulating in northeast China, and they pro- 
vide valuable epidemiological information for further studies of CCoV. However, further stud- 
ies are necessary to clarify the effect of co-infections and amino acid substitutions on the 
severity of the clinical symptoms of CCoV infections. 


Supporting Information 


S1 Fig. Nucleotide sequences of the partial M gene of CCoV strains identified in our study. 
(TXT) 


$2 Fig. Nucleotide sequences of the partial S gene of CCoV strains identified in our study. 
(TXT) 


S1 Table. Characteristics of the CCoV positive dogs, genotyping of CCoV strains, and 
amino acid substitution of the partial M protein in northeast China. 
(DOC) 


Author Contributions 


Conceived and designed the experiments: DBS LF. Performed the experiments: XYW (first 
author) CQL DHG SW YFG EYW ZHW XYW (fourth author). Analyzed the data: XWZ MJS 
QJL SYZ. Contributed reagents/materials/analysis tools: DBS. Wrote the paper: DBS. 


References 


1. Binn LN, Lazar EC, Keenan KP, Huxsoll DL, Marchwicki RH, Strano AJ. Recovery and characterization 
of a coronavirus from military dogs with diarrhea. Proc Annu Meet US Anim Health Assoc. 1974; 
78:359-366. PMID: 4377955 


2. Bandai C, Ishiguro S, Masuya N, Hohdatsu T, Mochizuki M. Canine coronavirus infections in Japan: 
virological and epidemiological aspects. J Vet Med Sci. 1999; 61:731-736. PMID: 10458093 


3. Decaro N, Buonavoglia C. An update on canine coronaviruses: viral evolution and pathobiology. Vet 
Microbiol. 2008; 132:221—234. doi: 10.1016/j.vetmic.2008.06.007 PMID: 18635322 


4. Ntafis V, Mari V, Decaro N, Papanastassopoulou M, Pardali D, Rallis TS, et al. Canine coronavirus, 
Greece. Molecular analysis and genetic diversity characterization. Infect Genet Evol. 2013; 16:129- 
136. doi: 10.1016/j.meegid.2013.01.014 PMID: 23410992 


5. Stavisky J, Pinchbeck G, Gaskell RM, Dawson S, German AJ, Radford AD. Cross sectional and longi- 
tudinal surveys of canine enteric coronavirus infection in kennelled dogs: a molecular marker for biose- 
curity. Infect Genet Evol. 2012; 12:1419-1426. doi: 10.1016/j.meegid.2012.04.010 22543007 PMID: 
22543007 


6. King AMQ, Adams Mu, Carstens EB, Lefkowitz EJ. Virus taxonomy: classification and nomenclature of 
viruses. Ninth report of the International Committee on Taxonomy of Viruses. Elsevier Academic 
Press, 2012, San Diego, CA. 


7. Le Poder S. Feline and canine coronaviruses: common genetic and pathobiological features. Adv Virol. 
2011; 2011:609465. doi: 10.1155/2011/609465 PMID: 22312347 


8. _Licitra BN, Duhamel GE, Whittaker GR. Canine enteric coronaviruses: emerging viral pathogens with 
distinct recombinant spike proteins. Viruses 2014; 6:3363-3376. doi: 10.3390/v6083363 PMID: 
25153347 


PLOS ONE | DOI:10.1371/journal.pone.0146975 January 15, 2016 9/11 


@-PLOS | ove 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 


Circulation of Canine Coronavirus in China 


Decaro N, Mari V, Campolo M, Lorusso A, Camero M, Elia G, et al. Recombinant canine coronaviruses 
related to transmissible gastroenteritis virus of Swine are circulating in dogs. J Virol. 2009; 83:1532— 
1537. doi: 10.1128/JVI.01937-08 PMID: 19036814 


Regan AD, Millet JK, Tse LP, Chillag Z, Rinaldi VD, Licitra BN, et al. Characterization of a recombinant 
canine coronavirus with a distinct receptor-binding (S1) domain. Virology 2012; 430: 90-99. doi: 10. 
1016/j.virol.2012.04.013 PMID: 22609354 


Costa EM, de Castro TX, Bottino Fde O, Garcia Rde C. Molecular characterization of canine coronavi- 
rus strains circulating in Brazil. Vet Microbiol. 2014; 168:8—-15. doi: 10.1016/j.vetmic.2013.10.002 
PMID: 24216489 


Decaro N, Cordonnier N, Demeter Z, Egberink H, Elia G, Grellet A, et al. European surveillance for pan- 
tropic canine coronavirus. J Clin Microbiol. 2013; 51:83-88. doi: 10.1128/JCM.02466-12 PMID: 
23100349 


Di Martino B, Di Felice E, Ceci C, Di Profio F, Marsilio F. Canine kobuviruses in diarrhoeic dogs in Italy. 
Vet Microbiol. 2013; 166:246—249. doi: 10.1016/j.vetmic.2013.05.007 PMID: 23806200 


Pinto LD, Barros IN, Budaszewski RF, Weber MN, Mata H, Antunes JR, et al. Characterization of pan- 
tropic canine coronavirus from Brazil. Vet J. 2014; 202:659-662. doi: 10.1016/j.tvjl.2014.09.006 PMID: 
25294661 


Terada Y, Matsui N, Noguchi K, Kuwata R, Shimoda H, Soma T, et al. Emergence of pathogenic coro- 
naviruses in cats by homologous recombination between feline and canine coronaviruses. PLoS One 
2014; 9:e106534. doi: 10.1371/journal.pone.0106534 PMID: 25180686 


Decaro N, Mari V, Elia G, Lanave G, Dowgier G, Colaianni ML, et al. Full-length genome analysis of 
canine coronavirus type I. Virus Res. 2015; 210:100-105. doi: 10.1016/j.virusres.2015.07.018 PMID: 
26221765 


Wang Y, Ma G, Lu C, Wen H. Detection of canine coronaviruses genotype | and Il in raised Canidae 
animals in China. Berl Munch Tierarztl Wochenschr. 2006; 119:35-39. PMID: 16450706 


Ma G, Wang Y, Lu C. Molecular characterization of the 9.36 kb C-terminal region of canine coronavirus 
1-71 strain. Virus Genes 2008; 36:491—497. doi: 10.1007/s11262-008-0214-4 PMID: 18306032 


Gao FS, Hu GX, Xia XZ, Gao YW, Bai YD, Zou XH. Isolation and identification of a canine coronavirus 
strain from giant pandas (Ailuropoda melanoleuca). J Vet Sci. 2009; 10:261-263. PMID: 19687628 


Li C, Wei S, Guo D, Wang Z, Geng Y, Wang E, et al. Prevalence and phylogenetic analysis of canine 
kobuviruses in diarrhoetic dogs in northeast China. J Vet Med Sci. 2015, Aug 7. [Epub ahead of print] 
PMID: 26256044. 


Geng Y, Guo D, LiC, Wang E, Wei S, Wang Z, et al. Co-Circulation of the Rare CPV-2c with Unique 
Gln370Arg Substitution, New CPV-2b with Unique Thr440Ala Substitution, and New CPV-2a with High 
Prevalence and Variation in Heilongjiang Province, Northeast China. PLoS One 2015; 10:e0137288. 
doi: 10.137 1/journal.pone.0137288 PMID: 26348721 


Pratelli A, Tempesta M, Greco G, Martella V, Buonavoglia C. Development of a nested PCR assay for 
the detection of canine coronavirus. J Virol Methods 1999; 80:11-15. PMID: 10403671 


Pratelli A, Decaro N, Tinelli A, Martella V, Elia G, Tempesta M, et al. Two genotypes of canine coronavi- 
rus simultaneously detected in fecal samples of dogs with diarrhea. J Clin Microbiol. 2004; 42:1797— 
1799. PMID: 15071054 


Erles K, Brownlie J. Sequence analysis of divergent canine coronavirus strains present in a UK dog 
population. Virus Res. 2009; 141:21-25. doi: 10.1016/j.virusres.2008.12.009 PMID: 19162099 


Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGAG6: Molecular evolutionary genetics analy- 
sis version 6.0. Mol Biol Evol. 2013; 30:2725-2729. doi: 10.1093/molbev/mst197 PMID: 24132122 


Buonavoglia D, Cavalli A, Pratelli A, Martella V, Greco G, Tempesta M, et al. Antigenic analysis of 
canine parvovirus strains isolated in Italy. New Microbiol. 2000; 23:93-96. PMID: 1094641 1 


Gémez MM, de Mendonga MC, Volotao Ede M, Tort LF, da Silva MF, Cristina J, et al. Rotavirus A geno- 
type P[4]G2: genetic diversity and reassortment events among strains circulating in Brazil between 
2005 and 2009. J Med Virol. 2011; 83:1093—-1106. doi: 10.1002/jmv.22071 PMID: 21503926 


Lau SK, Woo PC, Yeung HC, Teng JL, Wu Y, Bai R, et al. Identification and characterization of boca- 
viruses in cats and dogs reveals a novel feline bocavirus and a novel genetic group of canine bocavirus. 
J Gen Virol. 2012; 93:1573—1582. doi: 10.1099/vir.0.042531-0 PMID: 22495233 


Mesquita JR, Nascimento MS. Molecular epidemiology of canine norovirus in dogs from Portugal, 
2007-2011. BMC Vet Res. 2012; 8:107. doi: 10.1186/1746-6148-8-107 PMID: 22776749 


Choi S, Lim SI, Kim YK, Cho YY, Song JY, An DJ. Phylogenetic analysis of astrovirus and kobuvirus in 
Korean dogs. J Vet Med Sci. 2014; 76:1141-1145. PMID: 24784439 


PLOS ONE | DOI:10.1371/journal.pone.0146975 January 15, 2016 10/11 


®-PLOS | ove 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


Circulation of Canine Coronavirus in China 


Jeoung SY, Ann SY, Kim HT, Kim D. M gene analysis of canine coronavirus strains detected in Korea. 
J Vet Sci. 2014; 15:495-502. PMID: 25234323 


Soma T, Ohinata T, Ishii H, Takahashi T, Taharaguchi S, Hara M. Detection and genotyping of canine 
coronavirus RNA in diarrheic dogs in Japan. Res Vet Sci. 2011; 90:205—207. doi: 10.1016/j.rvsc.2010. 
05.027 PMID: 20557915 


Cavalli A, Desario C, Kusi I, Mari V, Lorusso E, Cirone F, et al. Detection and genetic characterization 
of canine parvovirus and canine coronavirus strains circulating in district of Tirana in Albania. J Vet 
Diagn Invest. 2014; 26:563—566. PMID: 24928599 


Decaro N, Desario C, Billi M, Mari V, Elia G, Cavalli A, et al. Western European epidemiological survey 
for parvovirus and coronavirus infections in dogs. Vet J. 2011; 187:195—199. doi: 10.1016/j.tvjl.2009. 
10.027 PMID: 19932978 


Pratelli A, Tinelli A, Decaro N, Cirone F, Elia G, Roperto S, et al. Efficacy of an inactivated canine coro- 
navirus vaccine in pups. New Microbiol. 2003; 26:151-155. PMID: 12737196 


Pratelli A, Tinelli A, Decaro N, Martella V, Camero M, Tempesta M, et al. Safety and efficacy of a modi- 
fied-live canine coronavirus vaccine in dogs. Vet Microbiol. 2004; 99:43—-49. PMID: 15019110 


Pratelli A, Martella V, Pistello M, Elia G, Decaro N, Buonavoglia D, et al. Identification of coronaviruses 
in dogs that segregate separately from the canine coronavirus genotype. J Virol Methods 2003; 
107:213-222. PMID: 12505636 


Decaro N, Martella V, Ricci D, Elia G, Desario C, Campolo M, et al. Genotype-specific fluorogenic RT- 
PCR assays for the detection and quantitation of canine coronavirus type | and type II RNA in faecal 
samples of dogs. J Virol Methods 2005; 130:72—78. PMID: 16024100 


Buonavoglia C, Decaro N, Martella V, Elia G, Campolo M, Desario C, et al. Canine coronavirus highly 
pathogenic for dogs. Emerg Infect Dis. 2006; 12:492—494. PMID: 16704791 


Decaro N, Buonavoglia C. Canine coronavirus: not only an enteric pathogen. Vet Clin North Am Small 
Anim Pract. 2011; 41:1121-1132. doi: 10.1016/j.cvsm.2011.07.005 PMID: 22041207 


Ntafis V, Mari V, Decaro N, Papanastassopoulou M, Papaioannou N, Mpaiziou R, et al. Isolation, tissue 
distribution and molecular characterization of two recombinant canine coronavirus strains. Vet Micro- 
biol. 2011; 151:238-244. doi: 10.1016/j.vetmic.2011.03.008 PMID: 21481551 


PLOS ONE | DOI:10.1371/journal.pone.0146975 January 15, 2016 11/11 


